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Abstract. There is growing evidence, based on both X-ray and radio observations of isolated neutron stars, that 
besides the large-scale (dipolar) magnetic field, which determines the pulsar spin-down behaviour, small-scale 
poloidal field components are present, which have surface strengths one to two orders of magnitude larger than 
the dipolar component. We argue in this paper that the Hall-effect can be an efficient process in producing such 
small-scale field structures just above the neutron star surface. It is shown that due to a Hall-drift induced 
instability, poloidal magnetic field structures can be generated from strong subsurface toroidal fields, which are 
the result of either a dynamo or a thermoelectric instability acting at early times of a neutron star's life. The 
geometrical structure of these small-scale surface anomalies of the magnetic field resembles that of some types of 
"star-spots". The magnetic field strength and the length-scales are comparable with values that can be derived 
from various observations. 
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1. Observational evidence 

There exists both observational and theoretical evidence 
of the existence of strong small-scale magnetic field struc- 
tures in the surface layers of isolated neutron stars. More 
and more, it becomes clear that neutron stars of nearly 
all ages must posses magnetic field structures, which are 
much more complicated than the simple assumption of a 
star centered magnetic dipole suggests. The polar surface 
magnetic field strength of such a dipole is conventionally 

estimated by Bd ~ dA-lO^^V Pp G, where standard neu- 
tron star quantities are assumed and the rotational period 
P (in seconds) and its temporal derivative P are taken 
from radio and/or X-ray timing observations. 

Recently, iBecker et al.l l)2003j) reported the (marginal) 
detection of an electron cyclotron line at about 3.3 keV in 
the Chandra-spectra of the millisecond pulsar B1821-24. 
The energetic location of that line corresponds to a surface 
magnetic field strength of 3 • lO^'^ G, while P and P yield 
a dipolar surface field strength of about 4.5 • 10^ G, almost 
two orders of magnitude lower. 

iHaberl et all l|2003 l) analysed the XMM-spectra of the 
isolated neutron star RBS 1223 and found evidence for 
a proton cyclotron absorption line in the energy range of 
0.1 . . .0.3 keV, which is consistent with a magnetic sur- 
face field strength of (2 ... 6) • lO^^ G. While P = 10.31 s 
seems to be a well settled value for the rotational period. 
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the determination of P is a much more delicate task in 
this case. Recent evaluations of XMM data indicate, how- 
ever, that the global dipolar field attributed to P and 
P is at least one order of magnitude smaller than the 
field corresponding to the proton cyclotron line (Schwope 
2003, private communication). Likewise, for the pulsar IE 
1207.4-5209 a surface field of - 1.5 • lO" G has been es- 
timated (Sanwal ct al. 2002), w hereas the estimate for its 
dipolar field is (2 ... 4) • lO^^ G ^Pavlov et al.ll2002(l . 

These X-ray observations indicate that, apart from a 
large-scale dipolar magnetic field which determines the 
spin-down behaviour of the pulsar, much stronger but 
short-ranged field components close to the neutron star 
surface exist. They do not affect significantly the braking 
of the star's rotation, but can affect the magnetospheric 
processes in the vicinity of the polar cap (defined by the 
totality of all open dipolar field lines). 

It is commonly accepted that pulsar radio emission 
is generated within a dense electron-positron plasma, 
the creation of which requires an ultra-strong poten- 
tial drop that accelerates charged particles along curved 
magnetic field lines. The observed phenomenon of drift- 
ing subpulses strongly suggests that this potential drop 
results from the deficiency in the actual charge den- 
sity wit h respect to the so-called co-rotational charge 
density tGoldreich fc JuHanI Il969l) just above the po- 
lar cap surface (see iGil eraLn200^ . The formation of 
suc h a charge-depleted reg i on (c alled 'polar gap' af- 
ter iRuderman k. Sutherlan who proposed it for 
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the first time) requires a strong and highly non-dipolar 
surface magnetic field, with radii of curvature much 
smaller than 10^ cm and magnitudes close to 10^'^ G, 
irrespective of the di polar field strength inferable from 
the pulsar spin-down i|Gil fc Mitrall200lt iGil fc Melikidzd 
l2002t iGil et ani2nn.'^ . Moreover, several periodicities ob- 
served in the phenomenon of drifting subpulses strongly 
suggest that the s ubpulse-associated plasnia fila ments 
(called 'sparks' after lRuderman fc Sutherland 1975') circu- 
late around a local magnetic pq le^ (;.Deshpande fc R,ankinl 
Il999[ l20nit [Gil fc Sendvkll20n.1) . Accordingly, the small- 
scale surface magnetic field anomalies are supposed to 
show spot-like structures allowing a persistent arrange- 
ment of drifting sparks (due to the well known E x B 
plasm a drift mechanism) in the form of qua si-annular pat- 
terns llGil fc Sendvkll20nflHFan et al.ll2nnih . 



2. The Hall-drift induced instability 

It is well known that the Hall-effect, probably via 
a cascade I Goldreich fc Reiseneggeil Il992j) . causes the 



generation of smaller scaled magnetic field components 
out of an existing large scale field. This process has 
been discussed for neutron s tar magnetic fields by a 
number of auth ors (see, e.g., 'Shalvbkov fc Urpinf 'l997l 
IVainsht ein et al.i 2000: HoUerbach & Riidiger 2002). Here, 
however, we want to refer to the fact that small- 
scale poloidal fields close to the neutron star surface 
can be generated from a subsurface toroidal magnetic 
field by a Hall-drift induced instability (HDI). The 
basic prerequisite for this instability is that a suffi- 
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strate, that strength and spatial structure of the HDI 
modes resulting from a realistic NS crust model are con- 
sistent with quantities derived from observations. 

The decay of a magnetic field in the almost crystallised 
crust of a NS is governed by 



B = curl ( Tj curl B + a(curl B x B)) 
div B = , 



(1) 



where the diffusivity rj and the Hall-parameter a deter- 
mine the intensity of Ohmic decay and Hall-drift, respec- 
tively, in their spatial distribution. Linearization of Eq. 
with respect to a background field Bq yields 

b = curl (77 curl 6) 

- curl (a (curl Bo X b -I- curl f) X Bo)) (2) 
div 6 

describing the evolution of a small magnetic perturba- 
tion b. For conditions as realised in the crusts of cool- 
ing isolated neutron stars, unstably growing perturbations 



^ A surface region with a magnetic flux prevailingly normal 
to the surface and with a certain degree of axisymmetry. 



b ~ b{r) exp 'yt , 7 > , can exist, with typical growth 
times 7^^ ^ lO'* "^ years. 

The existence of this instability has been shown in 
a slab geometry, approximating the neutron star crust 
geometry locally. Typical crustal density profiles were 
adopted. Boundary conditions are defined by assuming a 
perfect conductor being adjacent to the slab at the bot- 
tom and vacuum being adjacent at the top. That is, cor- 
responding to the superconductivity of the NS core, the 
magnetic and electric fields do not penetrate under the 
bottom of the slab whereas the magnetic field is contin- 
ued above the slab as a potential field approximating the 
conditions outside the NS. 

The background fields Bo we used are parallel to the 
slab; they vary strongly with depth and vanish at the bot- 
tom of the slab. Fields of that kind, especially toroidal 
ones, may exist as relics of a short, but very efficient con- 
yective dynamo phase acting in the proto-ne utron star 
ijThompson fc Duncanll9'93lUrDin fc Gii200.l() . However, 
it is perhaps more likely that such background fields are 
due to the action of a thermoelectric instability, which 
amplifies toroidal seed fields very e ffectively and sets on 
compu lsory when T^g/gsn ^ 100 l|Wiebicke fc GeppertI 
lll996() . Ts6 - surface temperature in 10^ K, gsi4 - sur- 
face gravity in 10^* cm s~^). Such conditions are met in 
almost all newly born neutron stars and are maintained 
up to an age of ^ 1000 years. Moreover, when a strong 
temperature gradient is established again by any external 
or internal process, the thermoelectric instability will be 
switched on again. It can thus, at later stages too, produce 
a toroidal background field which is capable of providing 
the conditions for the onset of the HDI. Hence, the initial 
strength of the background field is determined essentially 
either by the vigour of the dynamo or by the strength 
of the temperature gradient and may well exceed 10^^ G 
locally inside the crust. 

Of course, in any real situation the background field is 
not exactly parallel to the slab: Within the region we are 
interested in, the polar cap, a normal background field 
arises quite naturally. Its influence on the properties of 
the HDI can at the moment only be extrapolated from 
the special case of a homogeneous normal background field 
component in a homogeneous slab. Then, the tangential 
part of the background field rotates about a normal axis 
with an angular velocity proportional to the normal field 
component. The HDI continues to exist if the ratio of nor- 
mal to tangential background field components is not too 
large. If the normal field is small, what we want to assume 
here, the background configuration rotates slowly and we 
expect no major differences in the perturbation modes in 
comparison with those obtained for completely tangential 
background fields to appear. 

The HDI is described in detail for conditions realised 
in the crusts of isolated cooling neutron stars in RKG03 
and we refer the reader to that extended presentation. Its 
main result is that the HDI occurs for a variety of realistic 
crust conditions and that the obtained growth times are 
short enough to cause observable consequences. 
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The easiest way to get an idea on how the HDI acts 
is considering its analogy to field generation by hydrody- 
namic dynamo action. One major ingredient of dynamo 
models, e.g., those explaining the solar magnetic field, is 
shear motion (in axisymmetry typically occurring as dif- 
ferential rotation) which is capable of generating strong 
toroidal (in axisymmetry: azimuthal) fields from weak 
poloidal (in axisymmetry: meridional) fields by 'winding 
up'. In our context, the motion of the electron fluid, in 
which the field is partially frozen in, is well able to pro- 
vide shear (see the term a curlBo x 6 in Q). To get a 
successful dynamo, shear has to be completed by another 
effect generating in turn poloidal fields from toroidal ones. 
Mean-field dynamo theory has identified various possibil- 
ities for that, amongst which one (the so-called us x j 
effect) has exactly the same mathematical struct ure like 
the se cond Hall term acurlb x Bq in ((2Jl (see, e.g.. lRadieil 
Thus, a cycle of mutual amplification of poloidal 
and toroidal fields can establish resulting in a growing to- 
tal field. 



3. The small-scale poloidal field 

The structure of a typical^ unstable (i.e., exponentially 
growing) perturbation field (or eigenmode) is shown in 
Figs.nand|21(for detailed explanations see RKG03). While 
the background fields used in that paper are derived from 
the equatorial regions of global dipolar fields, we neverthe- 
less apply here the corresponding results in the vicinity of 
the polar cap region. There, a global poloidal field can 
no longer provide a suitable background field structure - 
a toroidal field is needed. Applicability of the results of 
RKG03 then requires that the toroidal background field 
near the pole and the poloidal one near the equator may 
have more or less similar radial profiles. (Note, that the 
properties of the HDI modes are only moderately sensi- 
tive with respect to the radial profile.) This seems to be at 
least possible in the case of convective-dynamo generated 
fields and even a reasonable assumption for thcrmoelectri- 
cally generated fields with not too low multipolarity. For a 
maximum strength of the toroidal background field in the 
order of about 10^* G the perturbation rises with a char- 
acteristic growth time of about 5 • 10"* years. We assume, 
that at saturation it may reach a significant fraction of 
the background field's strength. Then it is possible that 
its surface strength reaches 5 • 10^^ G and we scaled the 
perturbation field (the amplitude of which is not deter- 
mined by the linear Eq. (0)) simply to this value. 

The lifetime of such a perturbation is bounded from 
above by its ohmic decay time. Very likely, it is shorter be- 
cause the perturbation is subject to the Hall-effect which 
is supposed to accelerate the decay in general. With the 
conditions at the depth in which the perturbation currents 



^ Pandharipande-Smith equation of state, initial penetration 
density po 



10^'' g cm ^, cubic background field Bo with 



initial polar field strength 10 G, model age 3 • 10 years. 



circulate and their scales, we estimate the ohmic decay 
time of the field presented in Figs. ^ and |2] as > 10'' yrs. 

Both the major radial and tangential scales of the un- 
stable perturbations are given approximately by the radial 
extent of the background field. According to the assump- 
tions about the latter they are significant portions of the 
crust thickness (say < 50%) and scale with it. The crust 
thickness, in turn, is model dependent: Stiff equations of 
state (EOS) result in smaller compactnesses and hence 
larger star radii and crust thicknesses in comparison with 
softer ones. For example, the specific EOSs considered in 
RKG03, the stiff Pandharipande-Smith and medium soft 
Friedman-Pandharipande ones yield a crust thickness of 
3.8 km and 700 m, respectively. Generalising from these 
two cases we suggest that the scales of the eigenmodes are 
not very sensitive with respect to the star model. 

Let us stress that the spot-like structure presented in 
Figs, n^^nd 121 satisfies all conditions to create the vacuum 
gap and t o generate the elect ron -positron plasma within 
it. In fa,ct. lGi1 fc TviitT^ l|2nmh and lGil MefikidzfJ l|2nn2f> 
argued that the vacuum gap can only be formed under 
the so-called near-threshold conditions, when the surface 
magnetic field ^ O.IB, - 4.4 • lO^^ q addition, it 
requires small radii of curvature Rc "C 10^ cm. As can 
be seen from Fig. ^ in our case Rc < 10^ cm. In such 
a strong and curved surface magnetic field the magnetic 
pair creation via curvature and/or ICS (inverse Conipton 
scattering) photons is very efficient. 

Moreover, the magnetic field lines of the perturbation 
field converge at local poles (spot centers) as can be seen 
in Fig. 13 As already mentioned at the end of Sect. ^ 
this property of the surface magnetic field is crucial for 
the subpulse drift phenomenon. Of course, as an implicit 
assumption one has to keep in mind that just one spot 
structure plays the role of a local magnetic pole requiring 
that the canonical polar cap (i.e., the one formed by the 
dipolar field alone) coincides at least partially with this 
spot. The actual polar cap is defined by the open field 
lines of the superimposed global star-centered-dipole and 
local perturbation fields, along which an ultra-high accel- 
erating potential drop can exist, or, in other words, by 
those perturbation field lines which merge with the open 
dip olar field lines be fore entering the radio emission region 
fsee lGil et aLll2002l for details). 

When assuming that close to the surface the global 
dipolar field is at least five times smaller than the pertur- 
bation field the superposition of both fields would yield a 
field very similar to the one shown where only the symme- 
try between the regions with positive and negative normal 
flux is slightly disturbed. 

Note, that the spot-like structures of the surface field 
must not be conf used with the Deshpande-R ankin 'hot 
spots' or 'sparks ' jDeshpande fc RankinI [l999l) . The for- 
mer can be considered stationary within all the time scales 
relevant for observations and provide just the conditions 
for the development of sparks. By virtue of the E x B 
plasma drift the latter move around the local magnetic 
pole. As the appearance of the pure ExB drift is confined 
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Fig. 1. Structure of a typical magnetic field perturba- 
tion generated by the Hall-drift induced instability from a 
toroidal crustal background field having a maximum field 
strength of about 10^^ G. The growth time of the pertur- 
bation is of the order of 5 • 10'* years. The neutron star 
model behind relies on the Pandharipande-Smith equa- 
tion of state, its crustal thickness d ^ 3800 m. The vertical 
scale z is arbitrarily limited to z = d/\Q = 386 m, where 
the magnetic field strength is reduced by a factor of about 
3 compared with its surface value. The distance between 
adjacent spots in x-direction ~ 1.1 km, in y-direction 
~ 3 km (see Fig. |21 for that dimension). Colour encod- 
ing corresponds to the normal field component, < ■ >x 
denotes the average with respect to x. All magnetic field 
values are given in units of 10^^ G, Z4 is in units of 10^ cm 



to homogeneous fields only, additional drift components 
and/or acceleration prevent the charge carriers forming 
the sparks from moving exactly periodically along closed 
paths. So it remains to be examined whether our field 
spots are suited to reproduce the periodicities of the sub- 
pulse phenomenon in detail. 
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Fig. 2. The same magnetic field perturbation as in Fig.^ 
but viewed from above and zoomed into a smaller x in- 
terval. Note, that this field has to be superimposed on the 
global dipolar field having its pole in the vicinity of one 
of the spot centers to get the actual open magnetic field 
lines. 



Goldreich, P. & Reisenegger, A. 1992, ApJ, 395, 250 
Haberl, F., Schwope, A., Hambaryan, V., Hasinger, G., & 

Motch, C. 2003, A&A, 403, L19 
Hollerbach, R. & Riidiger, G. 2002, MNRAS, 337, 216 
Pavlov, G. G., Zavlin, V. E., Sanwal, D., & Triimper, J. 

2002, ApJ, 569, L95 
Radler, K.-H. 1969, Monatsber. Dt. Akad. Wiss., 11, 194 
Rheinhardt, M. & Geppert, U. 2002, Phys. Rev. Lett., 88, 

101103 

Rheinhardt, M., Konenkov, D., & Geppert, U. 2003, A&A, 

submitted, astro-ph/0310917 
Ruderman, M. & Sutherland, P. 1975, ApJ, 196, 51 
Sanwal, D., Pavlov, G. G., Zavhn, V. E., & Teter, M. A. 

2002, ApJ, 574, L61 



Geppert, Rheinhardt, Gil: Spots in Neutron Star Surface Magnetic Fields 



5 



Shalybkov, D. & Urpin, V. 1997, A&A, 321, 685 
Thompson, C. & Duncan, R. 1993, ApJ, 408, 194 
Urpin, V. & Gil, J. 2003, A&A, submitted 
Vainshtein, S., Chitre, S., & Olinto, A. 2000, Phys. Rev. E, 
61, 4422 

Wiebicke, H.-J. & Geppert, U. 1996, A&A, 309, 203 



